Interpretation of the angular dependence of the de Haas-van Alphen effect in MgB2 
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We present detailed results for the amplitude and field dependence of the de Haas-van Alphen 
(dHvA) signal arising from the electron-like tt sheet of Fermi surface in MgB2 . Our data and analysis 
show that the dip in dHvA amplitude when the field is close to the basal plane is caused by a beat 
between two very similar dHvA frequencies and not a spin-zero effect as previously assumed. Our 
results imply that the Stoner enhancement factors in MgB2 are small on both the a and n sheets. 



C3 

o 
o 

Oh 



I 

O 

o 



> 

in 
m 

O 



o 



X 
S3 



Soon after the discovery of superconductivity in 
MgB2, several groups reported calculations of its band- 
structure. The topology of the calculated Fermi surface 
was verified by measurement of the de Haas- van Alphen 
(dHvA) effect. The observed-^ dHvA frequencies were in 
very good agreement with the calculated^iii^ extremal ar- 
eas on each of the four sheets of Fermi surface. The mea- 
sured quasiparticle effective masses were also in excellent 
agreement with the calculated values, providing direct 
evidence for the predicted large difference in electron- 
phonon coupling constants between the a and tt bands. 

For fields less than 20 T, three dHvA frequencies have 
much larger amplitude than the others; two of these arise 
from the minimal and maximal extremal areas of the 
smaller of the two tubular a sheets, whereas the third 
originates from a neck on the electron- like tt sheet. In 
what follows we refer to these dHvA orbits as Fi , F2 and 
F3 respectively (a diagram showing the calculated Fermi 
surface and predicted dHvA orbits can be found in Rcf. 
i). 

InRef.H it was shown that the angular dependence of 
the dHvA amplitude for Fi and F2 could be adequately 
explained by the usual Lifshitz-Kosevich expression for 
the oscillatory torque Tqsc of a 3D Fermi liquid,^ but for 
-Fa it could not. There were two puzzling features. First, 
the amplitude of F3 showed a pronounced dip at ~ 76° 
(see Fig. which was attributed to a 'spin-zero' ef- 
fect. The Stoner enhancement deduced from the position 
of this dip is approximately two times larger than that 
predicted^ and four times larger than those measured for 
the a sheet orbits, Fi and F2. Second, we were unable 
to explain the angle dependence of the amplitude. In 
this paper we will present detailed measurements of the 
field dependent dHvA amplitude of F3 as a function of 
angle, and show that the feature previously ascribed to 
a 'spin-zero' is actually caused by a beat with another 
dHvA frequency. 

Quantum oscillations in the torque produced by a 
small single crystal of MgB2, (mass=5.6/ig) were mea- 
sured with a piezoresistive, doped silicon cantileverJ 
Changes in the resistance of the cantilever are directly 
proportional to the torque and were measured using an 
AC bridge techniqueiS The torque values are reported 



here in units of bridge resistance R, i.e, the off-balance 
voltage divided by the excitation current (the change in 
lever resistance is 4 times larger than this). We estimate 
that F = 10"i"R (F in Nm and R in Q.). The noise level 
is around 2mri or ~ 10~^^ Nm. The crystal of MgB2 
was grown by a high pressure synthesis routed and is the 
same as sample B in Refs.QandQ 

We interpret our data using the standard expression for 
the first harmonic of the oscillatory part of the torque for 
a 3D Fermi liquidSiia 
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where F is the dHvA frequency [F — (fi,/27re)A, A is 
the extremal orbit area in k-space]; A" = d^A/dk^ 
is the curvature factor and ip is the phase. Rd, Rt 
and Rs are the damping factors from impurity scatter- 
ing, temperature and spin splitting respectively. The 
Dingle factor, Rb = exp(— ^^^f ), where kp is the or- 
bitally averaged Fermi wavevectoriSiii and £ is the quasi- 
particle mean free path. The thermal damping factor, 
Rt = X/sinhA: where X = {2n'^k^m*T)/{heB), m* 
is the quasi-particle effective mass. The factor Rsc ac- 
counts for the additional damping when the sample en- 
ters the superconduct ing state, and was studied in de- 
tail in MgB2 m Ref. The spin splitting factor is 
given by Rs = cos[(7rgm_B[l + S])/{2me)] where 1 -I- S' 
is the orbitally averaged exchange-correlation (Stoner) 
enhancement factor, g is the electron g-factor, is the 
free-electron mass. When g{l + S){mB/me) equals an 
integer, Rs = and the spin- up and spin-down Fermi 
surfaces beat out of phase to produce a 'spin-zero' min- 
imum in the dHvA amplitude. If the location of any 
of these 'spin-zeros' can be measured, and mB{9,(j)) is 
known from band-structure calculations, then the Stoner 
factor on each orbit may be deduced. Note that Rg does 
not depend on B, so near to a 'spin-zero' the dHvA am- 
plitude is suppressed at all fields. 

Although it is possible to fit the data directly to Eq. 
it is often more illuminating to extract the field de- 
pendent amplitude A of the dHvA oscillations by fitting 
small sections of the data (comprising of 1.5 oscillations) 
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FIG. 1: Reduced dHvA amplitude {A) versus inverse field for 
several rotation angles. The lines are fits of the data to Eq. 

The data for 9 = 61.7° have been multiplied by 8 for 
clarity. 



to 

Vosc{B) = A sin + j + + 6 (2) 

(the linear term accounts for the slowly varying back- 
ground torque and magnetoresistance of the cantilever). 
We then divide A by B^Rt to give the reduced ampli- 
tude A, which in the absence of other effects (i.e., i?sc = 
l)H is proportional to the Dingle factor Rd- The quasi- 
particle effective mass m* in the expression for Rt was 
determined by measuring the temperature dependence of 
the dHvA amplitude [for F3, m* = (0.456 ± 0.005)TOe at 
e = 70.8°]. 

In Fig. ^ we show reduced dHvA amplitude A versus 
inverse field ('Dingle' plot) for orbit F3 at selected an- 
gles. For 9 = 61.7° and 71.4° the behavior is strictly 
exponential, A oc exp{—a/B). However, for 9 > 71.4°, 
a appears to increase rapidly and the Dingle curves be- 
come markedly non-exponential. The dHvA amplitude is 
strongly reduced, but unlike the expected behavior close 
to a 'spin-zero', the reduction is not uniform at all fields. 
For 9 < 72° the increase in the coefficient a with decreas- 
ing angle is given by a = (SO/sin^) x e{9). The sin 61 
factor arises from the cylindrical nature of this section of 
Fermi surface (the band mass and the dHvA frequency 
increase by the same factor). The factor e{9) was deter- 
mined experimentally by fitting the data for 9 < 72°; we 
find, e{9) ~ 1 -I- 3.7(1 — sin 6')^ reflecting an increase in 
scattering rate with decreasing 9. From this we estimate 
the quasiparticle mean free path (at 6* = 0°) on this orbit 
to be 660 A. 

We will show below that the most likely reason for the 
non-exponential Dingle curves is a beat with a second 



FIG. 2: Variation of the frequency difference AF and the 
relative amplitude /3 of the beat frequency as a function of d. 
The dashed line is Eq. I* = 0.12^ and AFo = 23.5 T. 



dHvA frequency. Before we present this analysis in detail 
we will briefly discuss two other possible explanations. 

The background and oscillatory torque cause deflec- 
tions in the cantilever so the measurements are not per- 
formed at strictly constant angle. This 'torque interac- 
tion' effecifi^ causes spurious generation of harmonics and 
mixing of frequencies, and because the amplitude is angle 
dependent it can also cause bending of the Dingle curves. 
We can quantify this effect by fltting the data to Eq. 
allowing for the change in dHvA frequency with fleld. To 
a good approximatioii F3 (9) — / sin{9) , and 9 is given 
by the measured angle plus the field dependent deflection 
of the cantilever, 9 — 9o + Xr{B) (A is a fitting constant 
which measures the 'stiffness' of the cantilever). For the 
data presented here we find that A ~ 0.04 degrees/fi;— 
The result of this analysis is that the torque interaction 
effect is much too small to explain the non-exponential 
behavior see in Fig. ^ 

Another possibility stems from a mosaic spread in the 
direction of the a or c axis of the crystal. Beats be- 
tween different dHvA frequencies originating from differ- 
ent part of the crystal can produce non-exponential Din- 
gle curves. -'^^ The size of the effect depends on the spread 
of the dHvA frequencies and in the present case would 
becomes largest as the field is rotated away from basal 
plane. Experimentally we find that curvature in the Din- 
gle plot only occurs for angles close to the basal plane, 
and therefore conclude that in our crystals the mosaic 
structure does not play a significant role. This conclu- 
sion is reinforced by the fact that the observed behavior 
is highly reproducible between different crystals. 

Simple trigonometric relations show that if two dHvA 
signals with frequencies F and F + AF and amplitudes 
in the ratio /3 (with identical Dingle factors) are added 
together, and the resultant fitted with a single frequency 
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F then the reduced amplitude A will vary like 
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In Fig. n we show that these equations provide a very 
good fit to the reduced dHvA amplitude data. In these 
fits there are three free parameters, Aq, AF and /3 (a was 
fixed at the values found by extrapolating from lower 9 
as described above) . The variations of AF and /3 with 9 
are shown in Fig. [3 The frequency difference increases 
approximately linearly with 9, AF ~ 1A{9 - 70°). The 
data are consistent with P being almost constant as a 
function of angle and then going rapidly to zero for 9 < 
73°. Note that as ^3(6*) ~ 2684/ sin(6'), the maximum 
value of AF corresponds to only ^ 0.4% difference in kp 
between the two orbits. 

A extremal dHvA orbit F3* with frequency close to F3 
was predicted from band structure calculations'^ and is 
shown in Fig. O Close to the two extremal orbits this 
sheet of Fermi surface is tubular, with only slight warp- 
ing. The difference in dHvA frequency between the ex- 
tremal orbits of a cylinder with simple cosine warpingiii^ 
is given by 
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Here /* is the fc-space distance between the minimum 
and maximum frequency orbits and Jo is the Bessel func- 
tion. For 9 > 65°, the detailed band structure determi- 
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F3 closely follows this equation 
with AFo ~ 57 T, and 1* ~ 0.13-^ (Ref. 19). For 

9 < 65°, AF3(6') = 0, i.e., only one extremal orbit is 
found. This reflects a difference between the actual Fermi 
surface topology and the cosine dispersion which leads to 



FIG. 4: Experimental dHvA torque amplitude versus 9 at 
fixed field (17.8 T < B < 18.0 T) for orbit F3. Data for two 
different in-plane rotation angles (j) are shown. Fits to the 
two frequency model [Eq. warped cylinder model [Eq. 

@] and single frequency model [Eq. iQ] are shown by the 
solid, dotted and dashed lines respectively. 



Eq. Q). The characteristics of this predicted second or- 
bit closely match our observations. The only significant 
difference is that the maximum frequency difference we 
find is around half the predicted value (i.e., the cylin- 
der is actually somewhat less warped than the calcula- 
tion). Experimentally, the Dingle curves are straight for 
9 < 72° (Fig.^, showing that the second orbit is absent 
for 9 ^ 72° in approximate agreement with the detailed 
calculation. 

We conclude that the dip in amplitude for orbit F3 is 
caused by a beat effect and not a spin-zero as previously 
assumed. From the data in Fig. 01 it seems likely that 
any spin zero would have to occur for 9 < 50° and hence 
[assuming for F3, iriB oc (sm9)~^] the Stoner enhance- 
ment factor on this orbit is less than 0.22. A summary of 
all the Stoner factors derived from dHvA measurements 
is shown in Table (the values for Fi and F2 are taken 
from a previous study^i on a different MgB2 crystal^"). 
The enhancement factors on both the a and tt sheets 
are somewhat smaller than calculations^ would suggest. 
By comparing measurements of the Pauli susceptibility, 
derived from conduction electron spin resonance experi- 
ments, to band-structure calculations of the total density 
of states at the Fermi level, Simon et alm^ found that 1+S 
(averaged over all Fermi surface sheets)was 0.86 ± 0.13, 
i.e., the average enhancement is small in agreement with 
our findings. 

Using the measured values of a, AF(6'), and f3{9) along 
with the calculated curvature factor A", we are able to 
calculate the expected angular dependence of the dHvA 
amplitude at fixed field using Eq. (QJ . Here we have fixed 
the constant of proportionality in Eq. Q to match the 
maximum in the experimental data, and we have also 
fixed 1 + = 1.2. As mentioned above, 1 -f S* must be 
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TABLE I; Summary of measured (1 + Scxp) and calculated- 
(1 + Scale) Stoner enhancement factors for the three orbits. 



Orbit 


Sheet 


1 ~t~ Scxp 


1 "h 5'calc 


Fi 


en 


1.07 


1.31 


F2 




1.12 


1.31 


F.i 


VTl 


< 1.22 


1.26 



< 1.22 and values less than 1.2 do not change the curve 
markedly. The agreement between the calculation and 
the data is good, although not perfect. Importantly, the 
calculation reproduces the key dip feature at ~ 76°. 

Actually, for such a flat portion of Fermi surface, the 
usual LK expression [Eq. |^] should be modified as the 
usual integral over the tube is no longer strongly peaked 
at the extremal orbits. For a simple cosine dispersion the 
curvature factor is replaced by^^ 
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A fit to the amplitude data with Eqs. Q and lO is 
shown as the dotted line in Fig.^ Here we have set I* = 
0.12-^ and AFq = 23.5 T to best fit both the AF{e) and 

^osc{0) data. The fit again reproduces the main features 
of the experimental data but is slightly worse than the 
discrete frequency fit. It can be seen (Fig. that the 



experimentally determined AF^^O) data does not quite 
follow Eq. Q , showing that the true local warping is not 
quite cosinusodial. The behavior below 9 = 74°, where 
the second frequency is not detected experimentally, is 
particularly badly described by this approximation. 

The dashed line in Fig. 01 shows the behavior expected 
without the beat effect but with 1 + 5 = 1.545, which is 
required to produce a 'spin-zero' dip at the same angle. 
Clearly the correspondence with the data is much worse, 
especially for angles close to the basal plane. 

For the series of sweeps reported in Figs. Q and 13 the 
in-plane angle (j> was estimated to be 8°. A second 
(less extensive) set of runs with ~ 0° shows very similar 
behavior. The main difference is that the dip occurs at 
lower 9 (see Fig. 0)). Repeating the above analysis on 
this set of data again shows a linear 9 dependence of AF 
but with a larger slope, dAF/d9 ~ 1.9 T/degree. If we 
use this value in our calculation we find the dip moves to 
lower 9, in accord with the experimental data. 

In summary, we have shown that the dip in the dHvA 
amplitude for fields aligned close to the basal plane of 
MgB2 is due to a beat between two very similar dHvA 
frequencies and not a spin-zero effect as previously as- 
sumed. The data imply that the Stoner enhancement 
factors on both the a and n sheets of Fermi surface are 
small. 

We thank J.R. Cooper for useful discussions and S. Lee 
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